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Nucleon Decay with Domain- Wall Fermions 
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'^RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 
''Physics Department, Columbia University, New York, NY 10027, USA 

We report on our on-going project to calculate the nucleon decay matrix elements with domain-wall fermions. 
Operator mixing is discussed employing a non-perturbative renormalization. Bare matrix elements of all the 
possible decay modes induced by the dimension-six operators are calculated with the direct method, which are 
compared with the indirect calculation using chiral perturbation theory. 



1. INTRODUCTION 

There are two methods to calculate nucleon de- 
cay matrix elemeirts using lattice QCD. One is 
the direct calculation of the matrix elements deal- 
ing with three- and various two-point functions 
|1I2IH| . The other is the indirect method where 
one relies on chiral perturbation theory whose 
low energy constants are calculated on the lat- 
tice [415111213] . The former method is preferable 
but requires huge computation, typically 10 times 
larger than the latter to achieve similar statistical 
accuracy. 

br a previous work ^ , we made use of good chi- 
ral symmetry of the domain- wall fermion with the 
DBW2 gauge actioir gj. The proton decay ma- 
trix element of the pion final state was calculated 
with both methods. The results are consistent 
with each other within the relatively large error. 
There we assumed the mixing of operators aris- 
ing from the explicit chiral symmetry breaking is 
negligible. Here we reconsider this issue by calcu- 
lating the mixing with a non-perturbative renor- 
malization (NPR). Also the matrix elements of 
the proton decay in case of the non-degenerate 
quark mass in the pseudoscalar state are calcu- 
lated allowing further discussion on the effective- 
ness of the chiral perturbation theory. 

All of our calculations are done on the 
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quenched DBW2 gauge configurations with size 
X Nr = 16^ X 32. The domain-wall height is 
M — 1.8, the fifth dimension size is Ls = 16 for 
the NPR and 12 for the matrix element. 

2. OPERATOR MIXING 

The nucleon decay is induced by the dimension- 
six operator which consists of one lepton and 
three quarks. The three quark operator, with 
which we calculate matrix element with initial 
nucleon and final pseudoscalar state, in general 
takes a form as 



(1) 



The superscript a specifies a combination of the 
gamma matrix (FF') listed in Tabled Operators 
are classified by the ordinary parity (V) and by 
the parity of switching (5) u and d quarks. The 
5 is a symmetry of the (lattice) QCD Lagrangian 
when quark masses are degenerate. The u, d, s 
do not necessarily mean the real flavor, but sim- 
ply stand for the different flavors. The other ar- 
rangements of flavors, such like O"^^ and O^us ^"^^ 
always expressed in terms of O^^^ by Fierz rear- 
rangement and/or by S. Operators with two fla- 
vors like Oudu follows the same argument. Lattice 
operators are renormalized as 
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For example, for operators in 7^+, sector, 
which will not mix the operators in the other sec- 
tors, a and h in eq. |21may be SS, PP, AA. The 



2 



Table 1 

Operator classification by parity and switching 
symmetry. 





s- 


s+ 




SS PP AA 


VV TT 




SP PS AV 


VA TT 



same 3x3 matrix Z applies to the V , S sector: 
SP, PS, -AV, too. 

For the NPR with RI-MOM scheme [TEl we 
calculate the Green's function of the operator 
with the three quark external states in the Lan- 
dau gauge, 

G''{xo,xi,X2,X3) = {0^j^^{xa)u{xi)d{x2)s{x3)).{3) 

We set xq = 0. Fourier transformation is per- 
formed on the other legs with the same momen- 
tum p, which are then amputated to obtain the 
vertex function, 



A^ip^) = F.T. {0,Xi,X2,X3)\ Amp. 



(4) 



Writing the tensor indices explicitly, the renor- 
malization condition of the RI-MOM scheme 
leads, 



(5) 



pa y— 3/2 ybc \ c ea6 

^ijk Pa &l' ^Hjk afl — " ' 

where Zq is the quark wave function renormal- 
ization i,j,k are color indices, a, (3 and 7, 6 are 
Dirac indices associated with F, F' respectively. 
The renormalization condition should be applied 
where the pion pole effect j8j is absent and the 
lattice artifact is small : Aqcd ^ \p\ ^ 
The projection matrix for example onto SS, 
is P^^ = e*J'=(C-i)'^"(5'5^/(6-4-4). We define the 
matrix M as, 

M'^'' — P.^jj^ ■ h'l^j, ^g, (6) 

which wiU give Zq^^iZ-^)''''. 

Fig. n shows the resulting at a quark mass 
TTif — 0.025 for the domain-wall fermion. 51 con- 
figurations have been analyzed . The renormal- 
ization factor of SS, SS and PP, PP are iden- 
tical within the error except for the vicinity of 
zero momentum, the others are zero within the 
two standard deviation level. Absence of the off- 
diagonal element and coincidence of the SS, SS 
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Figure 1. Elements of M for mj — 0.025 as a 
function of momentum squared. 



and PP, PP element are expected in the case of 
exact chiral symmetry. We have examined four 
different masses m/ = 0.025,0.04,0.055,0.07, 
where mass dependences are hardly seen in all 
range of momentum except near zero. 

As a result, the nucleon decay operator 

'^R/L-L = e'^Hu'^CPRfLd')PLS>' [S±P S-P] 

is renormalized multiplicatively without any mix- 
ing at practical level for this parameter of the 
domain-wall fermion. In the next section we 
show the results of the bare matrix elements with 
Lg = 12. The residual mass m^es = 1-6 MeV at 
Ls = 12 is still comparable to 0.7 MeV at Lg = 16 
[H] . Hence the effect of Lg on the renormalization 
factor is expected to be absent in this case, too, 
which needs to be checked by further performing 
simulation with Lc = 12. 



3. MATRIX ELEMENTS 

We follow the ref. |2] to calculate the nucleon 
decay matrix elements by the direct method. Our 
results for the decay to a meson with degener- 
ate quark mass are already shown in [3], where 
we have used 100 gauge configurations. We ex- 
tend our calculation to the non-degenerate quark 
mass case to treat the decay to kaon. Our masses 
are m/i,m/2 = 0.02,0,04,0.06,0.08 with m/i < 
mf2, where to/i is for the u, d quarks, m/2 for 
the u, d or s quark. The other parameters are the 
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-<itj(ud)„ujp> 
<Jt°|(ud)LuJp> 



<KJ(us)„uJp> 
<K°|(us)lUJp> 

-<K*|(us)„dJp> 
<K*|(us)LdJp> 

-<K*|(ud)„sJp> 
<K*|(ud)LSjp> 
<K°|(us)„dJn> 
<K°|(us)LdJn> 

<nl(ud)RULlP> 
<Ti|(ud)LUjp> 
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Figure 2. Summary Wo for both direct and indi- 
rect calculations with unrcnormalized operators. 

same as in Result of the relevant form factor 
Wq of the matrix element for various decay modes 
are shown in Fig. [3 We are assuming the SU(2) 
symmetry for the u and d quarks. There are other 
possible matrix elements, but they can be calcu- 
lated with the matrix elements in the figure when 
the SU(2) symmetry is intact. 

The indirect method uses the chiral pertur- 
bation theory, which relates the low energy pa- 
rameters in the chiral Lagrangian to the matrix 
elements. JLQCD extended the calculation by 
Claudson et. al JOl to all possible decay modes 
[2]. The low energy parameters are mass ratio 
of nucleon and baryon with strangeness S = —1, 
rriN/'mB, pion decay constant /, D and F pa- 
rameterizing the semileptonic decay of baryons, 
nucleon decay parameter a and /3. These param- 
eters are in principle calculable on the lattice. As 
we are using the chiral fermion, observables cal- 
culated on the lattice should obey the chiral per- 
turbation whose parameters are calculated on the 
lattice. However we use the experimental values 
except for a and The reasons are that we know 
the results of and gA C3 (= D + F) are 
anyway consistent with the experiment, and that 
we have no estimate of and individual D and 
F. Using the a and l3 in our previous study |3j 
and D = 0.8, F = 0.47 :i2, / = 0.131 GeV, 
mN/mB ~ 0.817, the matrix elements are cal- 
culated and shown in Fig. |21 The lattice cutoff 
calculated by p mass input, a^^ = 1.31 GeV is 



used to make dimensionless value. 

The indirect method gives quite good estimate 
of the matrix element. Consistency of the two 
• direct method within the error is in contrast to the re- 

° indirect suh of JLQCD 2^ with Wilson fermion. However, 

this could be caused by larger statistical error in 
our calculation. We need more statistics to be 
able to make more reliable statements. 

4. SUMMARY 

1 0.15 We have investigated the mixing of three- 

quark nucleon decay operators with domain-wall 
fermions by NPR. The mixing is practically neg- 
ligible at the parameter set similar to that where 
the matrix element are calculated. We have ex- 
tended our nucleon decay matrix element to those 
for the final kaon state matrix elements. The di- 
rect and indirect methods give consistent results 
within our statistical error. Performing NPR at 
the parameter of matrix element calculation to- 
gether with the matching calculation will give us 
the matrix element in MS. 
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